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Abstract 

Scattering of light by a random stack of dielectric layers represents a one-dimensional 
scattering problem, where the scattered field is a three-dimensional vector field. We 
investigate the dependence of the scattering properties (band gaps and Anderson local- 
ization) on the wavelength, strength of randomness and relative angle of the incident 
wave. There is a characteristic angular dependence of Anderson localization for wave- 
lengths close to the thickness of the layers. In particular, the localization length varies 
non-monotonously with the angle. In contrast to Anderson localization, absorptive 
layers do not have this characteristic angular dependence. 
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1 Introduction 

Scattering of waves in periodic structures (e.g. in a crystalline solid material) can be de- 
scribed by Bloch's Theorem [TJ [2], a theory that gives extended, propagating waves. A 
completely different situation appears if the periodicity of the scattering structure is dis- 
turbed by disorder: The scattered waves do not propagate any longer but become localized 
due to complex interference effects, provided that the disorder is sufficiently strong. This 
phenomenon, also known as Anderson localization, was originally proposed for quantum 
states [3] and later discussed in more detail in terms of a renormalization approach by Abra- 
hams et al. [4j. An important finding of the latter is that the tendency towards localization 
is stronger in low dimensions than it is in higher dimensions. According to this approach, 
quantum states are always localized in one and two dimensions, regardless of the strength 
of disorder. However, there are exceptions from this result. One case, where delocalized 
states can appear in the presence of random scattering in two dimensions, are relativistic 
(Dirac) states [5]. The main difference between nonrelativistic (Schrodinger) and relativistic 
(Dirac) states is that the former are scalar and the latter are spinor (i.e. vector) states, 
indicating that the dimensionality of the state (i.e. scalar vs. vector) plays a crucial role in 
the appearance of localization. 

Scattering of electromagnetic waves in a random ensemble of scatterers is a common 
problem in physics, biology, engineering, and astronomy. It can be used for a remote anal- 
ysis of complex objects like, for instance, the surface of a distant planet. A characteristic 
feature of light scattering is coherent backscattering (a comprehensive review can be found 
in [6]). In this context, a more exotic and, at least experimentally, less understood subject 
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is Anderson localization of light. This has been studied in terms of theory and experiment 
by a number of groups [TJ El EE [TOj HH [121 02] • m particular, the existence of Anderson 
localization for electromagnetic waves was discussed in Ref. [Hj. In contrast to the lo- 
calization of quantum states, the electromagnetic field is always a three-dimensional vector 
field. However, there are special situations in which the scattering process affects only one 
component of the electromagnetic vector field. An example is a wave scattered by a stack of 
layers with different refractive indices (cf. Fig. 1): If the wave vector of the incident wave is 
perpendicular to the layers, the vector components are scattered separately, and the result- 
ing scattering equation is a scalar one (the Helmholtz equation) which is formally equivalent 
to the Schrodinger equation of a quantum state. The Helmholtz equation has been widely 
used to study scattering in infinite and semi-infinite random media [151 EE El EE EH] • F° r 
the case of polarization, however, this can only be determined from the vector form of the 
Maxwell equation [20ll2l]. 

A particular case of a scattering medium is the layered system with a number of interest- 
ing features [2"2"l [231 EH ES] • The arrangements of layers is either periodic with an alternating 
refractive index or the refractive index changes randomly from layer to layer. Realizations of 
layered systems can be found in biological tissues, in the atmosphere, and in coated optical 
devices. The effect of Anderson localization then can be used to protect and cover an object 
if the latter is coated with random layers of non-absorptive materials. The advantage of this 
effect is that the electromagnetic waves are prevented to enter the coated object without 
absorbing the electromagnetic waves. 

An important aspect of the layered system is that the scattering is effectively a one- 
dimensional process. Since the vector components of the electromagnetic wave scatter inde- 
pendently only if the incident wave is perpendicular to the layers, by changing the direction 
of the incident wave one can couple the components of the electromagnetic wave in the scat- 
tering process. In other words, we can tune our scattering process from being a scalar one 
to being a vector one. Thus the scattering by layers will allow us to study the dependence 
of localization effects on the scalar or vector nature of the scattered field. 

In this paper we study the localization of electromagnetic waves by a stack of layers 
with randomly chosen refractive indices. The strength of localization is characterized by the 
Lyapunov exponent of the scattered wave in the direction perpendicular to the layers. It can 
be understood as the inverse of the localization length which can be measured in experiments. 
The advantage of considering the Lyapunov exponent in our calculations is the fact that it 
is believed to be a self-averaging quantity [26]. This enables us to perform the calculation 
without averaging over random ensembles. Our aim is to compare Anderson localization 
with absorption, where the latter is described by an imaginary part in the refractive index 
[27] . Both effects lead to an exponential decay of the intensity of light along the scattering 
process. The difference of their physical origin, however, should lead to different behavior 
with respect to wavelength and angle of the incident light. 



We consider a 3D system with random dielectric layers that are perpendicular to the x 
direction with refractive index n(x). For the electric field E(x,y,z) we can use the ansatz 
stationary plane wave in y and z direction: 



2 The Model 
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Using the short-hand notation E(x) = E'(x,k y ,k z ) and d x = d/dx, the electric field is 
determined by the Maxwell equation 



/ ky + k 2 ikyd x 
ikyd x k 2 d x 
\ ik z d x kyk z 



ik z d x 
kyk z 

k 2 y - d X 2 



E(x) = e(x)E(x) , e(x) = n{xf {u / cf . 



(2) 



For discrete layers and not too short wavelengths A = 2-kc/uj (i.e. for a wavelength A larger 
than the thickness of the layers), Eq. (j2j) can be written in a discrete form by integrating 
the Maxwell equation in x direction from a layer boundary x n to the next layer boundary 
x n+ \. Since E(x) does not change much within a layer of thickness d n = x n+ \ — x n if A > d n , 
we can replace its value for x n < x < x n+ i by 



E, 
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Moreover, the differential operator d x is replaced by a difference operator as 



d x E(x) -> — / d x E(x)dx 
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In other words, a; is replaced by the discrete coordinate n. For the following study we assume 
that the layers have the same thickness (but varying refractive index) and we choose this to 
be d n = d. Then all length scales are given in units of the thickness of the layers. Thus Eq. 
(121) can be written as 



^2(E n+ i + E n _i — 2E n ) + Ai(E n+1 — E n „x) + A E n = e n E n 
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This can also be expressed as a recursion relation of E n : 

C*n 1 ^4E n+1 + C n 1 BE n _i = E n 

with 

A = A 2 + A U B = A 2 -A 1 , C n = -A + 2A 2 + e n l. 
Multiplication of Eq. ([7]) by C n gives 

AE n+ i = C n E n — -BE n _i 

and a subsequent shift of n by 1 gives 

^4E n+ 2 — C n+ {E n+ i = —BE n . 

These two equations can be combined to a first order difference equation as 
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Now we multiply this equation from the left with 



A \ _1 _ / A' 1 
C n+1 A) ~ \A- l C n+l A- 1 A- 1 



/ n+l 

and introduce the new (six-dimensional) vector field 



(12) 



to write 

with the transfer matrix 

f 
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* n+2 = T n * n (14) 
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(15) 



The transfer matrix T n satisfies det(T n ) = 1. Eq. (fl4l) will be used as the starting point for 
our subsequent calculations. In particular, we can iterate this equation and obtain 

*2„+2 =f 2n ---f 2 f *o. (16) 

The special case of a plane wave propagating only in x direction (i.e. k y = k z = 0) leads 
to scalar equations for the electric field, since the Maxwell equation (j2J) decomposes into two 
independent scalar equations for E y and E z . For k y = k z = the transfer matrix becomes 
a 2 x 2 matrix, where the 3x3 matrices A and B are identical and equal to the scalar — 1. 
Moreover, C n also becomes a scalar: —2 + e n . As a consequence, the transfer matrix is 

Tn = ( h 1 . K = e n -2. (17) 



3 Results 

Elastic scattering of waves can be characterized by various physical quantities. A very useful 
quantity to characterize localization effects in our one-dimensional scattering geometry of 
the stacked layers is the Lyapunov exponent [26, 28J: 

7j = J™, -togfl-EnjI) U=x,y,z). (18) 

It measures the exponential decay of the magnitude of the wave due to scattering in the 
random medium. A non-localized wave has jj = 0, whereas jj 7^ describes a localized 
wave E n j. The larger |7j| the stronger the localization is. The self-averaging property of 
the Lyapunov exponent [26] is a crucial advantage for our numerical calculations: there is 
no need for averaging over an ensemble of random scatterers because this is achieved by 
choosing a sufficiently large stack of randomly chosen layers. 
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3.1 Analytic Results: Alternating Layers 



In the case of layers with alternating value of e n 

J e e = n 2 e uj 2 /c 2 for n even , , 

€n ~ \ e Q = n 2 u 2 /c 2 for n odd ' 1 U) 

the problem in Eq. (fT4l is translational invariant: T n = T. Therefore, it can be solved easily 
by diagonalizing the 6x6 matrix T. For k y = k z = (scalar case) we have only the two 
parameters h e ,h in the transfer matrix of Eq. ( fTTj) . Moreover, in Eq. ( TT6l) we need only 
the 2x2 version T n of Eq. ( TlTj) for even n: 

Thus the corresponding eigenvalues are 



A± = -1 + h a h e /2 ± ^-/i /i e + hlhl/A. (21) 

The eigenvalues are complex for < h Q h e < 4 with |A±| = 1. They represent propagating 
solutions 

f(X+) n \ T (e iRn \ T /nn . 

with the real "wave vector" k which satisfies 

cosk = — 1- (23) 
Solving this equation gives a dispersion oj 2 (k) (cf. Fig. 2) with a gap A, opening at K — ±7r: 



A 



1 1 



(24) 



For h Q h e > 4 we have pairs of real eigenvalues with A + A_ = 1 and no propagating 
solution. The six eigenvalues Xj (only \Xj\ is plotted) of the vector case (i.e. for k x = k y ^ 0) 
are shown in Figs. 3-6. Propagating solutions (i.e. |A 3 -| = 1) are found for different regimes 
oiu 2 /c 2 . 



3.2 Numerical Results 

In the case of scattering by random layers we rely on a numerical procedure. Performing 
such a calculation is relatively easy because we only need to multiply 6x6 (or 2 x 2 in the 
scalar case k y = k z = 0) randomly chosen transfer matrices. Then the Lyapunov exponent 
can be calculated from the product according to Eq. ffl8|) . However, when we perform many 
multiplications of transfer matrices, numerical accuracy plays a crucial role. The numerics 
is dominated by large eigenvalues but we are interested in eigenvalues Xj with \Xj\ ~ 1. The 
accuracy of the latter is suppressed by the large eigenvalues. In order to avoid this problem 
we apply a method which orthonormalizes the columns of the product matrix after a few 
multiplications via the Gram-Schmidt procedure. This process separates automatically the 
different exponentially growing contributions [29, 30J. Then the logarithms of the lengths 
of the vectors are stored. The Lyapunov exponents are given by the mean value of these 
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logarithms divided by the number of steps between orthonormalizations. Finally, the smallest 
Lyapunov exponent is stored. It was found for one-dimensional systems that the number of 
multiplications required for convergence is inversely proportional to the Lyapunov exponent 
and it is approximately given by [30] 

iVmax « 2(e 2 7 )-\ (25) 

where e is the relative accuracy. In order to reach convergency of e ~ 0.01 we calculate 
typically up to 10 6 layers (cf. Fig. 7). 

4 Discussion and Conclusions 

A stack of layers with alternating refractive indices n e and n Q presents an instructive example 
for the influence of the wave vector of the incident wave on the scattering properties. The 
relative angle y of the incident wave with the layers is given through 

cosy = y/l - c 2 {k 2 + k 2 z )/u 2 . (26) 

If we start with the scalar or y = 0) we find two bands of propagating 

waves with a gap at the boundaries of the Brillouin zone (cf. Fig. [2]). As soon as we 
introduce a small nonzero wave vector parallel to the layers (i.e. k y = k z ^ or cosy < 1), 
the transfer matrix becomes six- dimensional. The gap between the two bands of propagating 
solutions remains almost unaffected (cf. Fig. [3]). Besides the two eigenvalues of the scalar 
case, related to the propagating solution, there are also eigenvalues \Xj\ ^ 1 which are not 
related to solutions of our scattering problem. 

For y > the band of propagating solutions (i.e. for eigenvalues \Xj\ = 1) persists, 
together with the band gap (cf. Fig. [5] for k y = k z = 2) but the band gap disappears for 
n Q = n e (cf. Fig. Hj), as we anticipate from the scalar case. However, in contrast to the scalar 
case, there is only a tiny band gap for n ^ n e , as shown in Fig. [6J 

In the case of random layers we consider randomly independent fluctuations of the re- 
fractive index in Eq. (J2]) 

n 2 n = 1 + ah n , (27) 

where h n is a random variable, distributed uniformly on the interval [—0.5, 0.5]. The param- 
eter a controls the strength of the random fluctuations. We find that (Anderson) localization 
is increasingly efficient with decreasing wavelength A = 2ttc/uj = 2ix/k. This is plausible, 
since a wave with large wavelength experiences the randomly fluctuating refractive index as 
an average refractive index. In Fig. [8] the Lyapunov exponent ^ x is plotted as a function 
of a, where 7 X increases monotonously with a and k. On the other hand, the behavior of 
the Lyapunov exponent as a function of the angle y is not so obvious. For uj/c ~ 1 it has 
a non- monotonous behavior (cf. Fig. [9]): the Lyapunov exponent decreases linearly with 
cosy, reaches a minimum near cosy = 0.8 and then increases. For uj/c ~ 1.2 there is well- 
developed minimum on the interval 0.6 < cosy < 0.8 , as shown in Fig. [TUJ For larger values 
of u/c (i.e. shorter wavelengths) the Lyapunov exponent is almost constant as a function 
of y. The non-monotonous behavior of j x with respect to the angle y might be related to 
the fact that the mechanism of localization is due to interference of different parts of the 
scattered wave. The distance between two sucessive scattering events, which happen at the 
interfaces of the layers, depends of the angle y. Changing y implies a change of the distance 
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between the scattering events and, therefore, allows us to change between constructive and 
destructive interference. 

Absorption is another mechanism which leads to an exponential decay of the wave E n j. 
It can be included in our Maxwell equation (j2J) by adding an imaginary part to the refractive 
index [27]: 

nl + iv- (28) 

In a homogeneous medium absorption should not be very sensitive to a change of the angle ip. 
Assuming identical absorptive layers (i.e. a = 0) we find, in contrast to the case of random 
layers, no characteristic variation of the Lyapunov exponent (cf. Fig. [TT1) . This effect 
allows us to distinguish Anderson localization and absorption in a scattering experiment by 
measuring the variation of the Lyapunov exponent as a function of the angle (p. 

In conclusion, we have studied the propagation of an electric field in a stack of equally 
thick layers. Layers with alternating refractive index do not allow propagation for some 
wavelengths by opening a band gap. Layers with randomly changing refractive index, on the 
other hand, have no propagating solution but experience Anderson localization. The strength 
of the latter depends crucially on the wavelength and the relative angle of the incident wave. 
In contrast to Anderson localization, absorption does not show a characteristic dependence 
on this angle. 

References 



[i] 

[2] 
[3 
[4 

[5 
[7 



Bloch F. Z Physik 1928:52:555-600. 

Kittel C. Quantum Theory of Solids. New York: J. Wiley, 1987. 
Anderson PW. Phys Rev 1958:109:1492-1505. 

Abrahams E, Anderson PW, Licciardello DC, Ramakrishnan TV. Phys Rev Lett 
1979:42:673-676. 

Ziegler K. Phys Rev Lett 1998:80:3113-3116. 

Mishchenko MI, Travis LD, Lacis AA. Multiple Scattering of Light by Particles: Cam- 
bridge: Cambridge University Press, 2006. 



John S. Phys Rev Lett 1987:58:2486-2489. 

De Raedt H, Lagendijk A, de Vries P. Phys Rev Lett 1989:62:47-50. 
[9] Wiersma DS, Bartolini P, Lagendijk A, Righini R. Nature 1997:390:671-673. 
[10] Scheffold F, Lenke R, Tweer R, Maret G. Nature 1999:398:206-207. 
[11] Wiersma DS, Rivas JG, Bartolini P, Lagendijk A, Righini R. Nature 1999:398:207. 
[12] Ziegler K. J Quant Spec Rad Transf 2003:79-80:1189-1198. 
[13] Milner V, Genack AZ. Phys Rev Lett 2005:94:073901-1-4. 
[14] Figotin A, Klein A. J Opt Soc Am A 1998:15:1423-1435. 



7 



[15] Akkermans E, Wolf PE, Maynard R. Phys Rev Lett 1986:56:1471-1474. 

[16] van Tiggelen BA, and Tip A. J Phys I France 1991:1:1145-1154 

[17] Klyatskin VI, and Saichev AI. Uspechi Phys Nauk 1992:162:161-194 

[18] Freilikher V, Pustilnik M, Yurkevich I. Phys Rev Lett 1994:73:810-813. 

[19] Chang SH, Cao H, and Ho ST. IEEE Journal of quantum electronics 2003:39:364-374. 

[20] Stephen MJ, Cwilich G. Phys Rev B 1986:34:7564-7572. 

[21] Ozrin VD. Phys Lett A 1992:162:341-345. 

[22] Zhang Z-Q. Phys Rev 1995:B52:7960-7964. 

[23] Feng Y, Ueda K. Optics Express 2004:12:3307-3311. 

[24] Hu L, Schmidt A, Narayanaswamy A, Chen G. Journal of Heat Transfer 2004:126:786- 
792. 

[25] Bertolotti J, Stefano G, Wiersma DS, Ghulinyan M, Pavesi L. Phys Rev Lett 
2005:94:113903-1-4. 

[26] Deych LI, Zaslavsky D, and Lisyansky AA. Phys Rev Lett 1998:81:5390-5393. 

[27] van de Hulst HC. Light Scattering by Small Particles. New York: Dover Publications, 
1981. 

[28] Lifshitz IM, Gredeskul SA, and Pastur LA. Introduction to the Theory of Disordered 
Systems. New York: J. Wiley, 1988. 

[29] Pichard JL, and Sarma G. J. Phys. C: Solid State Physics 1981:14:L127-L132. 

[30] MacKinnon A, and Kramer B. Z Physik 1983:B53:1-13. 



8 



Figure 1: Scattering of an electromagnetic wave on layers with refractive index e n . 



CM 
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Figure 2: Dispersion (w/c) 2 (/t) of alternating layers with refractive indices n e = 1.000, n Q = 
1.025 and k y = k z = (scalar case). There is a gap opening at the edges of the Brillouin 
zone. 
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Figure 3: Eigenvalues of the 6x6 transfer matrix for an alternating array of layers with 
the refractive indices of Fig. 2 and k y = k z = 0.001 (close to the scalar case). The gap is 
not affected by the very small change of k y = k z in comparison with Fig. 2. 
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Figure 5: The same as Fig. 4 but for n e = 1.000, n = 1.005 and k y = k z = 2. A propagating 
solution has a gap between uj 2 9.9c 2 and 10.0c 2 (see inset). 
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Figure 6: The same as Fig. 5 but for n e = 1.000, n = 1.025. A new gap for the band of the 
propagating solution develops near oo 2 xs 0.85c 2 . 
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Figure 7: Scattering by random layers: Lyapunov exponent for E x as a function of the 



number of random layers N for k z = k y 
random variable with —0.5 < h„ < 0.5. 
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Figure 8: Behavior of the Lyapunov exponents for various wave vectors k = k y = k z and 
e n — e o(l + othn). h n is the random variable of the previous figure. e = 2 (for k = 0), e = 5 
(for k = 1.5), eo = 8.5 (for k = 2), eo = 13 (for k = 2.5). The wavelength in the direction 
perpendicular to the layers is 27r//c, in units of the thickness of the layers. 
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Figure 9: Behavior of the Lyapunov exponents for random layers with a = 1 as a function 
of the angle ip of the incident wave and u 2 /c 2 . 




Figure 10: Behavior of the Lyapunov exponents for a system of random layers a — 1. 
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Figure 11: Behavior of the Lyapunov exponents for absorbing layers with 77 = 0.1 
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